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ABSTRACT 

The restframe UV-to-optical flux ratio, characterizing the "UV upturn" phenomenon, is potentially the most 
sensitive tracer of age in elliptical galaxies; models predict that it may change by orders of magnitude over the 
course of a few Gyr. In order to trace the evolution of the UV upturn as a function of redshift, we have used the 
far-UV camera on the Space Telescope Imaging Spectrograph to image the galaxy cluster CL0016+16 at z = 0.55. 
Our 25" x 25" field includes four bright elliptical galaxies, spectroscopically confirmed to be passively evolving 
cluster members. The weak UV emission from the galaxies in our image demonstrates that the UV upturn is 
weaker at a lookback time ^5.6 Gyr earlier than our own, as compared to measurements of the UV upturn in 
cluster E and SO galaxies at z = and z = 0.375. These images are the first with sufficient depth to demonstrate 
the fading of the UV upturn expected at moderate redshifts. We discuss these observations and the implications 
for the formation history of galaxies. 

Subject headings: galaxies: evolution — galaxies: stellar content — ultraviolet: galaxies 



1. INTRODUCTION 

Because they are composed of old, passively evolving pop- 
ulations, elliptical galaxies offer great promise for tracing the 
evolution of the Universe. Whether elliptical galaxies formed 
through hierarchical merging or monolithic collapse, one of the 
major goals in extragalactic studies is the determination of the 
"redshift of formation," zf , that marks the age where most of 
the stars in early-type galaxies formed. Recent studies of galaxy 
clusters out to z ~ 1 indicate that most of the star formation had 
to be completed at high redshift (z ;> 3), followed by quiescent 
evolution thereafter (Stanford, Eisenhardt, & Dickinson 1998; 
Kodama et al. 1998). 

The UV upturn is a sharp rise in the spectra of E and SO 
galaxies shortward of restframe 2500 A. It provides a sensitive 
tracer of age for the oldest stars in these galaxies, and can poten- 
tially constrain zf ■ Traditionally characterized by the 1550 — V 
color, the UV upturn in local galaxies originates in a population 
of hot horizontal branch (HB) stars and their U V-bright progeny 
(see Brown et al. 2000 and references therein). As first noted 
by Greggio & Renzini (1990), the UV upturn should evolve 
rapidly with age, through the evolution of HB morphology and 
the main sequence turnoff mass. Although all models of ellipti- 
cal galaxy evolution predict a rapid evolution in the UV upturn 
(e.g., Tantalo et al. 1996), the timing for the UV upturn onset 
depends strongly upon model parameters. 

We have been undertaking a series of observations to trace 
the evolution of the UV upturn as a function of redshift. Faint 
Object Camera (FOC) observations of Abell 370 (z = 0.375) 
provided the first detection of far-UV emission from quiescent 
elliptical galaxies above z = 0.05 (Brown et al. 1998). The very 
strong UV emission found at z = 0.375 suggests no evolution 
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in the UV upturn between our own epoch and one 4 Gyr earlier, 
a finding inconsistent with some models of galaxy evolution, 
and apparently consistent only for high values of zf ■ Here, we 
describe observations that trace the UV upturn to higher red- 
shift. We used the far-UV camera on the Space Telescope Imag- 
ing Spectrograph (STIS) to measure the UV emission from four 
giant elliptical galaxies in the cluster CL0016+16 at z = 0.55. 
Ground-based spectroscopy confirms their passive evolution 
and cluster membership (Dressier & Gunn 1992); morpholog- 
ical classification comes from Wide Field Planetary Camera 
2 (WFPC2) imaging in the F555W and F814W bands (Smail 
et al. 1997). We assume the currently popular cosmology of 
SIm = 0.3, f^A = 0.7, and// e = 67 km s Mpc~', although our 
results are more sensitive to zf than to the assumed cosmology. 

2. observations 

Using the STIS far-UV camera on 19 Aug and 21 
Aug 1999, we observed a 25" x 25" field centered at 
RA(J2000)=0 /, 18 m 33.2 J and Dec(J2000)=16°26'9.7" in the 
cluster CL0016+16. The total exposure was 27892 sec; 10 
frames were taken in two visits, with dithering by 6 pixels to 
smooth out small-scale detector variations. We used the crystal 
quartz filter (F25QTZ) for a bandpass that spans 1450-2000 A, 
thus reducing the sky background from geocoronal O I and 
Lyman-a to negligible levels at little cost in galaxy light, as the 
redshift (z = 0.55) puts the Lyman limit at the short wavelength 
cutoff. Because the long wavelength cutoff of this bandpass is 
due to detector sensitivity, red leak is also negligible; we avoid 
the problematic red leak and red grating scatter that hampered 
earlier attempts to measure the UV upturn at moderate redshift 
(e.g., Windhorst et al. 1994). The photometric calibration is re- 
liable at the 0.15 mag level (Baum et al. 1998). Recent efforts 
have revised the calibration slightly; we assume the latest revi- 
sion, with an accuracy that should be better than 0.15 mag. For 
reference, we assume that a flat spectrum of 1 .036 x 10~ 16 erg 
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Fig. 1- A schematic of the STIS image. Each bright cluster E galaxy 
is denoted by a cross, photometric aperture, and sky annulus. When 
there is a strong dark background, its shape follows the contours shown 
(contour levels are at 2x, 5x, and 10 x the lowest dark rate). A non- 
cluster spiral (asterisk) was used for astrometry and registration; it is 
shown in the inset at 4 x magnification, as it appears in the STIS (left) 
and WFPC2 (right) images. Triangles denote ~ 8a detections, and 
diamonds mark other objects from the Smail et al. (1997) catalog at 
/8i4 < 24 mag. Galaxies appearing in the Dressier & Gunn (1992) cat- 
alog are labeled, along with their morphology (Smail et al. 1997). 

sec -1 cm -2 A -1 produces one count per sec in the STIS band- 
pass. A full description of the instrument can be found in 
Woodgate et al. (1998) and Kimble et al. (1998). 

Unfortunately, the mean temperature of the far-UV camera 
has increased since STIS commissioning, and the dark back- 
ground increases as a function of temperature. This increase ap- 
pears as a "glow" that is strongest in the upper left-hand quad- 
rant of the detector (Figure 1), where the dark rate can be 20 
times higher than the nominal 6xl0 -6 cts sec -1 pix~'. The 
glow varies from weak to strong in our 10 frames. 

The Galactic foreground extinction toward our field was 
thought to be E(B — V) = 0.022 mag at the time these ob- 
servations were planned, based upon H I maps (Burstein & 
Heiles 1984). The subsequent release of a new extinction map 
(Schlegel, Finkbeiner, & Davis 1998), based upon an IRAS 
map of the Galactic dust, revises the extinction estimate to 
0.057 mag. The increased dark background and higher extinc- 
tion both work to decrease the sensitivity of our observations 
to far-UV emission, but do not preclude constraints on the UV 
upturn. When accounting for this extinction, we assume the 
empirical model of Cardelli, Clayton, & Mathis (1989), which 
has been parameterized and tested from the UV to IR. 

Our STIS field includes four giant E galaxies (see Figure 1); 
three are detected at >2.5er significance. We detect three 
other objects from the Smail et al. (1997) catalog at a signif- 
icance of ~ 8cr. We use the bright background spiral, DG244 
(z = 0.6577; Dressier & Gunn 1992), to confirm the registra- 
tion of our frames, and to determine the positions of the fainter 
galaxies. The accuracy of relative astrometry within WFPC2 or 
STIS images is better than 0.1", but the accuracy of absolute 
astrometry can be worse than 1", so the positions of objects in 



the STIS field are determined relative to this bright spiral. Most 
of the objects in our chosen field are too red and faint to be 
detected by STIS; their non-detection does not suggest an un- 
expected lack of sensitivity, and all detected and non-detected 
objects are consistent with the expected range of spectral en- 
ergy distributions (SEDs). STIS repeatedly images UV-bright 
stars (e.g., globular cluster NGC6681), and no drop in sensitiv- 
ity is apparent in observations before and after our observation 
of the cluster CL0016+16. 

3. DATA REDUCTION 

The STIS data was reduced via methods nearly identical to 
those used for the Hubble Deep Field South, to which we re- 
fer the interested reader for more details (Gardner et al. 1999). 
In brief, the images were processed via the standard pipeline, 
excluding the dark subtraction, low-frequency flat field correc- 
tion, and geometric correction. Dark frames from July through 
August of 1999 were processed in the same manner, and those 
with a strong glow (> 2 x 10~ 5 cts sec -1 pix~') were summed 
and fit with a cubic spline to produce a profile appropriate for 
our CL0016+16 observations (the shape of the glow changes 
slowly with time, and so recent darks are required for this fit). 
A flat component and glow component to the dark background 
were then subtracted from each CL0016+16 frame, and then the 
flat field correction was applied. The frames were registered by 
integer shifts and summed via the DRIZZLE package. The pix- 
els in each frame were weighted by the ratio of the exposure 
time squared to the dark count variance, including a hot pixel 
mask. The algorithm weights the exposures by the square of 
the signal-to-noise (S/N) ratio for sources that are fainter than 
the background, thus optimizing the summation to account for 
the temporal and spatial variations in the dark background. The 
statistical errors in the final drizzled image (cts pix~'), for ob- 
jects below the background, are given by the square root of the 
final drizzled weights map scaled by the exposure time. Geo- 
metric correction, which tends to smear out the pixels through 
non-integer shifts, was not applied. For small dithers and large 
extraction apertures, it is not needed for object registration. 

WFPC2 F8 14W data for the same field were obtained from 
the STScI archive and reduced via standard techniques, includ- 
ing cosmic-ray rejection and masking of problematic pixels. 
The total exposure in the F814W image is 16800 sec. The 
summed image was used to determine the restframe optical- 
band flux for each elliptical galaxy, and to determine the rela- 
tive positions of objects in the STIS field. 

4. PHOTOMETRY 

We performed aperture photometry on the far-UV frames us- 
ing IDL, taking a weighted average of the flux within a 16-pixel 
(0.4") radius, and a weighted average of the flux within a sky 
annulus of radii 80 and 100 pixels. The aperture includes the 
bright core of each galaxy, as they appear in the WFPC2 data. 
The aperture is small enough to minimize the background in the 
far-UV measurement, but large enough to avoid significant er- 
rors in the expected encircled energy (due to the uncertainty in 
the position of the galaxies); it is also relevant to measurements 
made of local galaxies, which use a nuclear aperture (see §||). 
The weighting used the map of statistical errors (see §[|), in the 
sense that pixels were weighted less if they had less exposure 
(due to masked pixels) or higher dark count rates. This weight- 
ing did not weigh by counts in the data frame, which would 
obviously bias the photometry toward pixels with more source 
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Table 1: Photometry 
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a Smail et al. (1997). 
b Dressier & Gunn (1992). 

counts. We determined the local sky value from the mean in- 
stead of the median in the sky annulus, because most pixels are 
ones and zeros. Note that alternate methods for determining the 
sky background (e.g., fitting a surface to the local residual back- 
ground) yield results well within the 1 a photometric errors. The 
positions of objects in the STIS image were determined from 
the geometrically-corrected WFPC2 F8 14W frame, relative to 
the spiral galaxy DG244 (see Figure 1). Calculation of position 
in the STIS frame involves a rescaling, rotation, translation, and 
geometric distortion, using the geometric distortion coefficients 
of Malumuth (1997). We tested our positioning algorithm us- 
ing WFPC2 and STIS images of the globular cluster NGC6681, 
which was observed with STIS numerous times before and af- 
ter our CL0016+16 observations at various roll angles and po- 
sitions. Absolute positional accuracy for objects in the STIS 
frame is 1-2 STIS pixels, well within the 16-pixel radius used 
for aperture photometry. 

For the WFPC2 F814W frames, we performed aperture pho- 
tometry with the IRAF package PHOT, using a 4-pixel radius 
(0.4"), and a sky annulus of radii 20 and 25 pixels, thus match- 
ing the photometry done in the STIS image. Note that this 
aperture size would produce encircled energy agreement at the 
5% level for point sources, and better agreement for extended 
sources (Robinson 1997; Holtzman et al. 1995), so the un- 
certainty in encircled energy contributes less than 0. 1 mag to 
1550 — V . We do not perform an aperture correction because 
we are only interested in colors, not the absolute fluxes. 

Photometry was performed on all objects with 7gi4 < 24 mag 
in the Smail et al. (1997) catalog, plus a faint object (#696) that 
is obvious in the STIS image but faint in the WFPC2 frames. 
Table 1 gives the photometry for the four giant elliptical galax- 
ies, plus three objects detected at ~ 8<r in the STIS image. 

5. INTERPRETATION 

The UV upturn is traditionally characterized by the restframe 
1550 — V color (see Burstein et al. 1988). Conversion to rest- 
frame 1550 — V from our observed bandpasses depends upon 
the assumption of an SED. Brown et al. (1998) used the spectra 
of three local elliptical galaxies (NGC1399, M60, and M49) to 
convert their observed FOC bandpasses to restframe 1550 — V, 
and we used the same spectra here. We redshifted the spectra 
of NGC1399, M60, and M49 to z = 0.55, and then applied a 
foreground reddening of E(B — V) = 0.057 mag (Schlegel et 
al. 1998). We then used the IRAF CALCPHOT routine to cal- 
culate the STIS and WFPC2 countrates, giving: 

1550 — V = —2.5 x \og]o(cps FUV /dpsfsuw) — 3.8 mag, 



where cpsfuv is the STIS countrate (cts sec -1 ), and dpsf^w 
is the WFPC2 countrate (DN sec -1 ). This conversion gives the 
restframe 1550 — V values shown in Table 1. We also calculate 
the 1550 — V assuming the 3cr upper limit to the far-UV flux. 
Note that the 1550 — V values in Table 1 would be redder if the 
foreground extinction is closer to E(B — V) = 0.022 mag (see 
§||), or if the far-UV SEDs are dominated by post-AGB stars 
(see below). The 1550 — V color for the sum of the giant ellip- 
tical galaxy photometry gives an average measurement of the 
UV upturn in our sample. 

Although we are looking to a significantly earlier epoch, the 
assumption of a model SED over an empirical SED makes lit- 
tle difference for the restframe V, because the observed F814W 
band overlaps with restframe V. For example, assuming that 
E galaxies are ^5 Gyr old at z = 0.55, we calculated the 
FEl AW observed — V rest frame color using the Bruzual & Chariot 
(1993) instantaneous burst SED of age 5 Gyr, and the NGC1399 
SED, assuming z = 0.55 with a foreground reddening of E{B — 
V) = 0.057. The difference in F814W oisen , e£ ; - V restframe for the 
two SEDs was only 0.02 mag. The assumed SED has somewhat 
more effect on the far-UV. Our STIS bandpass spans restframe 
935-1290 A, and the 1550 — V color is based on the average 
flux from 1250-1850 A. However, we note that no model SED 
has been tested in the far-UV for passively-evolving galaxies at 
moderate redshift, and thus we prefer empirical SEDs. Younger 
galaxies should have a larger contribution from relatively hot 
post-AGB stars, instead of HB stars (see Brown et al. 1997), so 
our conversion might systematically produce a bluer 1550 — V 
than the true restframe 1550 — V. Alternatively, if younger 
galaxies are dominated by even cooler HB stars than those in 
local galaxies, we are calculating a redder 1550 — V than the 
true value. 

Another source of systematic error is the aperture size. Lo- 
cal E galaxies have been measured through 14" diameter aper- 
ture (1.5 kpc at the distance of Virgo), sampling only the nu- 
clear flux (Burstein et al. 1988). The z = 0.375 measurements 
were made through a 1.82" diameter aperture (9.8 kpc at the 
distance of Abell 370), including all of the galactic light de- 
tected by the FOC (Brown et al. 1998). We have measured the 
flux from our z = 0.55 galaxies through a 0.8" diameter aper- 
ture (5.4 kpc at the distance of CL0016+16), enclosing the light 
from the galaxy cores. Our aperture is 3.6 times larger than 
that used for local galaxies, but 1.8 times smaller than that used 
for the Abell 370 galaxies. Ohl et al. (1998) found that local E 
galaxies usually (but not always) become redder in 1550 — V at 
increasing radius; the color gradient is not large (0.1-0.5 mag) 
over radii 1-25", and the color of the integrated light in an in- 
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creasing aperture will change even less, because most of the 
light comes from the central ~ 2 kpc even in our larger aper- 
tures. Thus, the aperture effects should be small, but perhaps 
non-negligible (~ 0.1 mag). 

6. DISCUSSION 

The average restframe 1550 — V color for the four giant E 
galaxies at z — 0.55 is much redder than that observed in clus- 
ters at z — and z = 0.375. At the 3a upper limit to the far-UV 
flux, none of the z = 0.55 galaxies are as blue as the strong 
UV-upturn galaxies observed locally (e.g., M60 or NGC1399). 
Ideally, one would want to compare galaxies at similar veloc- 
ity dispersion (cr v ) in each epoch, because <r v can be measured 
in a model-independent manner, and because the UV upturn 
strongly correlates with a v locally (Burstein et al. 1988). Mea- 
surements of <7,, are unavailable for some of the z = 0.375 galax- 
ies and all of the z = 0.55 galaxies, but these galaxies were se- 
lected from the brightest and largest in each cluster, and would 
likely show strong UV upturns if observed locally. The faint 
far-UV emission from the z = 0.55 galaxies thus demonstrates 
that the UV upturn is sensitive to age. We plot these UV upturn 
measurements in Figure 2, along with the models of Tantalo 
et al. (1996), which trace the chemical evolution of elliptical 
galaxies under the assumption of gas infall. Note that the scat- 
ter in local UV upturn measurements is affected by the wide 
range of sizes and luminosities in the sample; the moderate- 
redshift samples are more homogeneous, but the measurements 
themselves are less statistically significant. 

In the Tantalo et al. (1996) models, the UV upturn appears at 
an age of ^6.5 Gyr, becomes quite blue by 9 Gyr, and is flat 
thereafter. Assuming the models are correct and that the oldest 
stars in cluster E galaxies all formed at a common redshift, our 
measurements imply zf ~ 4. If elliptical galaxies form through 
hierarchical merging, the age of the oldest stars may predate 
such merging; alternatively, the age of the oldest stars may re- 
flect the age of the galaxies themselves, if they were assem- 
bled through monolithic collapse. However, the Tantalo et al. 
(1996) models rely on several parameters (e.g., time of onset 
for galactic winds, efficiency of the star formation rate, accre- 
tion timescale, etc.) that can be tuned to delay or accelerate the 
onset of the UV upturn. In this sense, our results say more about 
the evolution of HB morphology than zf ; if galaxies form their 
oldest stars at z > 4, our data imply that a blue HB population 
cannot arise until ages greater than 7 Gyr. 

Although the elliptical galaxies in our program are chosen 
through ground-based spectroscopy that confirms their passive 



evolution, the strong UV upturn found at z = 0.375 could the- 
oretically be due to residual star formation instead of evolved 
populations. This star formation would have to cease by the 
present epoch, because the UV emission from local cluster el- 
liptical galaxies is clearly due to HB stars instead of star for- 
mation. If all clusters are alike, our new z = 0.55 observations 
reinforce the case against star formation at z = 0.375, because 
such star formation would likely be even stronger at a larger 
lookback time. Measurements at z = 0, 0.375, and 0.55 demon- 
strate that the UV upturn rises sharply over a few Gyr, and then 
levels off. This line of inquiry into the evolution of galaxies is 
still in the early stages, and further observations are needed over 
a range of redshifts. Such measurements should become much 
more feasible in HST Cycle 10, once a cooler is connected to 
STIS (reducing the dark background to minimal levels). 

z 
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Fig. 2- Evolution of the restframe 1550 — V color according to the 
models of Tantalo et al. (1996) (shaded), viewed as a function of look- 
back time, assuming a reasonable set of cosmological parameters and 
three different epochs of galaxy formation (labeled). The spread in 
1550 - V is bounded by models at M = 3 x 10 12 and M = 10 12 M Q . 
Note the sudden onset of UV emission caused by the appearance of 
metal-rich hot HB stars. The points represent the mean 1550 — V color 
measured for quiescent E and SO galaxies in clusters at z — (Virgo, 
Fornax, and Coma; Burstein et al. 1988), z = 0.375 (Abell 370; Brown 
et al. 1998), and z = 0.55 (CL0016+16). The error bars at z = and 
z = 0.375 give the rms in galaxy colors (photometric errors are neg- 
ligible); the error bar at z — 0.55 gives the photometric error for the 
summed galaxies (rms in the galaxy colors is smaller). 
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